Tunneling in quantum wires I: Exact solution of the spin isotropic case. 
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We show that the problem of impurity tunneling in a Luttinger liquid of electrons with spin is 
solvable in the spin isotropic case (g a = 2, g p arbitrary). The resulting integrable model is similar to 
on ■ a two channel anisotropic Kondo model, but with the impurity spin in a "cyclic representation" of the 

on : quantum algebra su(2) q associated with the anisotropy. Using exact, non perturbative techniques 

we study the RG flow, and compute the DC conductance. As expected from the analysis of Kane 
and Fisher we find that the IR fixed point corresponds to two separate leads. We also prove an 
exact duality between the UV and IR expansions of the current at vanishing temperature. 
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I. INTRODUCTION. 



| Quantum impurity problems have attracted constant attention recently. The reason is, that the underlying physics 
!h ■ is highly non trivial, that the models are manageable technically despite the presence of very strong interactions, 
and that there are many practical/experimental applications. The latter include the Kondo effect, quantum dots, 
dissipative quantum mechanics, impurity tunneling in fractional quantum Hall devices and in one dimensional quantum 
wires. 

The physical variables are usually the temperature, the external field, the bulk interactions, the nature of the 
impurity (which can either be a charge or spin impurity), and the way it is coupled to the bulk degrees of freedom. 
A standard approach to computing physical quantities in these models would be to use perturbation theory in the 
bulk-impurity coupling for example. This perturbative expansion is most of the time done around the UV fixed point, 
and usually leads to limited results. Since all the problems of interest are essentially one dimensional (or can be 
reformulated as such), more powerful methods can however be used. For instance, in recent years, a conformal field 
theory description has allowed a much better understanding of the fixed points and their vicinity [jjj . 

Another method, which is a priori less general, is the use of quantum integrability, pioneered in the study of the 
■ Kondo problem 0] . Integrability is a powerful, non perturbative tool, which is not limited to the vicinity of the fixed 
points but allows a complete description of the properties all the way along the renormalisation group flow, from 
the UV to the IR fixed points. Technical difficulties have for a long time restricted the use of integrability to the 
computation of thermodynamic, equilibrium properties. However, recent progress, in particular in the understanding 
of form factors (^] and integrable description of conformal field theories ||,^| , have also made possible the computation 
of dynamical, out of equilibrium properties for many quantum impurity problems |||7] . 

A still open challenge in that field is the problem of a single (charge) impurity in a one dimensional quantum wire, 
where the electrons are described by a Luttinger liquid with charge and spin degrees of freedom. This problem was 
studied using the renormalisation group and perturbation theory by Kane and Fisher and Furusaki and Nagaosa 
I ||. There, a complete picture of the flows generated by the RG was obtained in terms of the "<?" factors, g p and g a , 
describing the interaction of bulk charge and spin degrees of freedom. The existence of tentalizing new fixed points, 
which would be partly transmitting and partly reflecting, was in particular conjectured. Steps were taken to identify 
O ■ these fixed points using conformal field theory |Io[ , but the problem remains open. 

In this paper, we are addressing this impurity problem from the point of view of integrability. The case of an 
impurity in a spinless electron gas could be mapped on the boundary sine-gordon model ||, which is a completely 
integrable quantum field theory [ p~Tj | , formally analogous to the anisotropic Kondo problem, but with the impurity 
spin in a cyclic representation of the quantum algebra su(2) q Jl2| . Here, we show that the case of electrons with spin 
maps, when g a = 2, onto a model analogous to the anisotropic two channel Kondo model, again with the impurity 
spin chosen in a cyclic representation of the quantum algebra su(2) q . This allows us to compute the relevant physical 
properties. 

II. THE MODEL, INTEGRABILITY, AND THE S MATRIX. 

In order to fix notations, let us repeat some of the basic definitions in Q. In the one dimensional Luttinger liquid, 
bosonisation of the fermionic operators is accomplished via 
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with /i =|, j. The fields M and M have commutation relations 

[^(x),9^(x')] = id fi , fl >9(x-x') (2) 

from which we can devise two different representations of the Luttinger liquid (our conventions follow Q). Here, we 
will work in the so called O-representation. Changing basis to the charge and spin degrees of freedom 



we obtain the action 



e p = e T + e i , e CT = e T -e i; (3) 

S = S p + S a (4) 
= f dxdy^-[(d x Q p ) 2 + (d y Q p ) 2 } + -L[(d x 9„) 2 + (d y G a ) 2 }. 

where y is the imaginary time, the spin and charge velocities have been normalised such that v a = v p = 1. In this 
convention, the g-factors have value g a — g p = 2 for a non-interacting system. The case of quantum wires with SU (2) 
symmetry corresponds to g a = 2. 

The electric and magnetic conductances of this system follow directly from Kubo's formula 

e 2 

Gp/a = -fr9p/<7- (5) 

As discussed in |D|, this formula might not always be physically relevant, because of the coupling to the reservoirs. It 
has been proposed that G p = 1%- might always hold, whatever the interactions in the quantum wire. Of course, this 
does not mean that the scattering through the impurity will be trivial in that case; it will indeed, still be controlled by 
the action written above, the only difference being in the meaning of the "applied voltage" . Here, we will simply follow 
Q for uniformity of notations. Physical results for different reservoirs configurations can be recovered by rescaling 
the voltage. 

In the presence of a charge impurity at the origin x — 0, the hamiltonian gets an additional piece 

5H = J dxV(x)(^ + ^i). (6) 

Here V(x) is a potential which has essentially zero measure outside x = 0. Under the bosonisation rules sketched at 
the beginning, this leads to the change of action 

SS = A J dy cos ^/nQ p cos \fixQ a , (7) 

where A ~ V(2kp ). In the physical case, there is no apparent symmetry that allows us to restrict to the action ^ - 
the most general form of the perturbation is 

5S = Jdy J2 ^^e^ n ? e ? +n ^ (8) 

where the A's are real couplings. The renormalisation group equations read at first order 

d*>n p ,n a ( n 2 nl \ 

- S -=^- T fl "-T fo J A "- B -- (9) 

As mentioned before, in this paper we restrict to the case g a = 2, which is the one relevant for physical quantum wires 
with unbroken SU(2) spin symmetry. Then, beside the perturbation with n p = n a = 1, the operators n p = 2,n ff = 
is relevant for g p < 1. The operator with n a = 2, n p — is always irrelevant. We will thus further restrict in the 
following to the (probably not physically relevant) case where only An 7^ or the case where g p > 1, so we have to 
take into account the term An only, which we simply call A in what follows. 
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In H, the equivalent problem for spinless electrons was exactly solved. The solution required a folding to transform 
the impurity problem into a boundary problem, and then used recent results on boundary integrable quantum field 
theories jll| together with the massless scattering approach. The same folding can easily be accomplished in the 
problem with spin. First, it is convenient to rescale the fields, writing the action as 

S = S p + S a + \ J dy cos y/Tfg^Qp cos V2irQ a (10) 

s„ = J dxd y ^[(d x ep) 2 + (d y e^) 2 }. (ii) 

We then introduce odd and even fields 

% = ^ [%,l(x, y) + %A-*> y)} (12) 

% = -j=[Q», L (x,y)-G l >, R (-x,y)]. 

With this, the interaction at x = involves only the even fields which have Neumann boundary conditions. The odd 
fields, having Dirichlet boundary conditions, completely decouple and do not interact. Also, the even field, as defined 
above, is left moving and we can "fold" using 

e fl , L = e e (x + y), x<o 

e^ R = e e (-x,y), x<0. (13) 
We can now express everything in terms of the fields 9 which is defined on the negative axis, the action becomes 
S = S + B 

= f dx I dy^J2 [(W 2 + (W 2 } + A f rfycos J^^(O)cos0F^(O). (14) 

The next step in J(| was to use the integrability of the corresponding boundary quantum field theory - the boundary 
sine-Gordon model (BSG). In the present case with spin, the problem, involving two fields, is a priori more complicated. 
Indeed it is unlikely that the problem is soluble for general values of g pi g a . However, in the case g a = 2, the problem 
can be solved, as we now demonstrate. 

The trick to solve a boundary problem such as the one under consideration here, is to introduce a more general 
integrable problem with both bulk and boundary perturbations. The bulk term is chosen such that the theory remains 
integrable thus defining a basis of quasiparticles bulk excitations. These excitations are interacting through an elastic, 
two particles S matrix derived from the constraints of unitarity, crossing-symmetry and the Yang-Baxter equation 
jl4]]. Since the theory with boundary is also integrable, there are also strong constraints on the boundary reflection 
matrix R leading to a very simple scattering of these quasiparticles at the boundary. When one lets the bulk coupling 
go to zero, these become massless quasiparticles which scatter in a simple way at the boundary. Taking the massless 
limit is by now a standard procedure and has been discussed in J^|,[^| . 

In our case, we need to find the appropriate bulk perturbation for the problem. The first natural idea to handle 
(Q) is to introduce as a bulk perturbation 

Sb = A / dx I dy cos 2\[k9o cos \p2mg~pQp (15) 

J — oo J 

Note the doubling of the cosine's arguments compared with (|l4|): this is because, quite generally, an integrable 
perturbation of the form f(0 L )f(0 R ) corresponds to an integrable perturbation f(0 R ) at the boundary. Here the 
function / is an exponential (to make the action real, appropriate combinations of exponentials were taken), and, 
using Neumann boundary conditions, we wrote cos#(0) = cos2# fl (0). The problem (15]) is however not integrable for 
generic values of g p (integrable sub-varieties have been identified in [p7|^9[ and will be the subject of a sequel to this 
paper), so this approach is not satisfactory. Fortunately, we can introduce another bulk perturbation, 

S B =k( dx I dy cos2y^9^ cos2y^6 R cos v / 2TTg p 9 p (16) 
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This corresponds to the function / being a cosine function. This problem is actually well known to be integrable. To 
see this, we can refermionize the 9 a part. The term cos2y / 7r0^ being the sum of a an exponential and its conjugate 
becomes now a Majorana fcrmion, so we have, equivalently, 



— oo 



S B = A / dx dy Wcos v /2ng p p (17) 



while the free part can be written as a free bosonic hamiltonian, plus the sum of two decoupled free Majorana fermions, 
only one of them appearing into the perturbation 

So = J dx J dy \[{d x e p f + (d v p ) 2 } + - VdV + x Bx - xdx- (18) 

Factoring out the other fermion (x), we end up with a standard integrable c = 3/2 analog of the sine-Gordon model 
|p3|. The boundary perturbation is also integrable. In terms of fermions, it can be written as 



A J dy cos y^f^,(0) [*(0) + *(0)] a(y) + \ J dy [*t(0) + ad] . (19) 

The last part simply enforces appropriate UV boundary conditions . The quantum operator a anticommutes with the 
fermions, and satisfies a 2 = 1 (see pCl| ] for a detailed discussion of bosonization in the presence of a boundary). 
In the bulk, the theory can be shown to be integrable by exhibiting a set of non-local conserved currents 

j± = ijfe ±ie p/V^, J± = -qle^p/V^ (20) 

These currents in turn define conserved charges with su(2)~ commutation relations. They commute with the charges 
constructed from the following local currents 

G = ^d6 p , G = l>d6 p (21) 

which are also conserved by the perturbation. Requiring that these symmetries are realized on the asymptotic states 
of the theory leads to a bulk S matrix factorized into a a su(2)~, sine-Gordon type part describing the scattering of 
a doublet, S u d, and a supersymmetric part, which is basically a kink S matrix with 3 vacua [ ^0[ . 

It must be mentioned here that there are some subtle issues regarding the integrability of the previous field theory - 
for simplicity we discuss only the bulk case here, the boundary case being quite similar. Although naive power counting 
following the arguments of |^l| suggests that the theory is integrable at any order, the fact that the fermions have 
integer dimension makes the problem more complicated. It is generally believed that the integrable theory requires an 
additional term of the form A 2 cos 2^/27rg p 0, as can be checked in the semiclassical case p2|] , or by imposing N = 1 
supersymmetry. The universal properties of the model ( |l7j ) are the same as the one with this additional coupling. 
Indeed, consider for instance the computation of the free energy perturbatively. If a is the short distance cut-off, its 
expansion involves the coupling A(R/a)2~ 29p for every insertion of the first operator, and the coupling A 2 (i?/a) 1 ~ 89p 
for every insertion of the second. For g p positive, the integrals are divergent. They could either be rendered finite by 
dimensional regularization (akin to taking g p negative), or by using the cut-off. In the latter case, the divergent part 
will only contribute to the regular, bulk part (see eg p3fl). Writing the free energy per unit length / = f reg — Tf, 
letting T = X/R — > and A — ► such that x = AR 1 ~ 29p remains finite, we obtain a function f(x) determined fully 
by insertions of the first operator, not of the second one. 

In fact one can avoid the appeal to this discussion of the bulk integrable theory, and N = 1 supersymmetry, by 
addressing the problem from the two channel Kondo model point of view. It is well established that the anisotropic 
Kondo problem is integrable for any number of channels and impurity spin; the case at hand corresponds simply to a 
particular impurity spin - a cyclic representation of the quantum algebra su(2) q . In particular, for the anisotropic two 
channel Kondo model, the S-matrix of the excitations in the bulk is precisely the one mentionned above (S u d® Skink)- 
This identification in addition gives the impurity S-matrix immediately - this is discussed in more details in appendix 
B. 

Recall that the solution of the massless problem proceeds by first describing the bulk problem in terms of massless 
relativistic excitations - the limit A — > of the excitations of (p7|). The basic massless particles are then either right 
or left moving, with dispersion relation e — p (e — —p), and they have a spin quantum number and a kink quantum 
number. We denote them by the doublet (u 0)0 ±i, d a>a ±i). The label a denotes one of the 3 possible vacua; a kink 
connects adjacent vacua, and is thus designated by a pair of labels (see figure 1). We use for the kinematic properties 
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of these basic particles the same parametrization, e = me , where m has the dimension of a mass, although the theory 
is massless. 




FIG. 1. Structure of the vacuas for the RSOS sector. 

The particles have factorized scattering, with trivial LR scattering. The LL and RR scattering are described by 
the bulk S matrix with tensor product structure 

S = S u d <g> Skink (22) 

Here S u d is the usual soliton-antisoliton sine-Gordon S-matrix, corresponding to the coupling ^ = — j^, ie a quantum 
group symmetry with deformation parameter q = — e~ J7r / 7 , where ^e. = 2 ^ +2 ^ , so the total dimension of the boundary 
perturbation (using that g a — 2) is d = | + 2 (y+2) = t+5- ^he S-matrix is the well known Restricted Solid on 
Solid S-matrix for a model with three vacua. There may also be bound states, depending on the value of 7. There 
are no such bound states if 7 is larger than one. When 7 is smaller than one, there are n of them, where n is the 
integer part of 1/7. 

The integrability of the model with boundary interaction translates into the fact that massless particles are reflected 
with no particle production, in a way described by a reflection matrix R, solution of the boundary Yang-Baxter 
equation. By analogy with the higher spin two-channel Kondo model (see appendix B), we find 

R = Ru,d ® 1 (23) 

where R u ,d is the reflection matrix of the boundary sine-Gordon problem at the foregoing coupling [3. It is enough 
for our purpose to recall the physical amplitudes 

Kl 2 = — J {9 _ lnTB y \<\ 2 = 1 - Kl 2 (24) 

I + 6 7 

Here Tb is an energy scale related with the coupling A by Tb oc X^ 1 /^ 1 ) oc A 7+2 . Physical quantities to be discussed 
below will be expressed in terms of the squares of R-matrix elements, and therefore expand in powers of T^ 1 . On 
the other hand, these quantities expand in powers of the square of the amplitude in front of the leading irrelevant 
operator determining the approach to the IR fixed point. If the IR fixed point is approached along fiO, with O of 
dimension d, then a has dimension 1 — d. Hence we find 2(d — I) = -orrf=- + I = = h + ^rp^ = — + — . 
Following the discussion in M , this indicates that the IR fixed point corresponds to "disconnected leads" , as expected 
in that domain of the parameters. 

III. DC TRANSPORT PROPERTIES 

The whole logic of |6) can now be implemented. To compute DC properties, it is enough to treat the massless 
particles as free ones from the point of view of impurity scattering, ie use a Landauer Biittiker type formula. The 
massless particles are quantized using TBA equations. The system of equations, for 7 an integer, is, setting / pj — e ej 
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e 3 (9, V/T) = 8^ -K*J2 N i* ln \ l + e^e-oWW] . (25) 



where the kernel K(9) = ^ cosh 9) » anc ^ -^j'fe ^ s the incidence matrix of the following diagram (which is found by 
"gluing" the RSOS and the non-diagonal sine-Gordon TBA diagrams) 



FIG. 2. TBA diagram at zero temperature. 

i.e. Njk = 1 if the nodes j and k are connected, and otherwise (in particular Njj = 0). The chemical potential 
vanishes, except for the end nodes, where p = ±( 7 + 1)^, V the applied voltage. Here, all the nodes but the one 
with j = 2 correspond to pseudoparticles, which are necessary to diagonalize the S-matrix. 

A. Closed form solution at T = 

The simplest case is when T = 0. Then, a Fermi sea of antisolitons is formed when a voltage is applied. All the 
scattering involves these antisolitons only, so the results can simply be derived for any value of 7. Introduce the 
quantity e, solution of 

<o) J A <£(e - e'W)de' = \-^e 6 (26) 

Here, $ is a kernel with Fourier transform 



1 1 sinh(7 - 1)ttw/2 

4 cosh 2 iru j 2 2 cosh ttuj/2 sinh jttuj/2 



$ H = ; + ; — w (27) 



It is convenient to write it in the form 

* = : Bi^( 7 + 2W2 (28) 

4 cosh itu) j 2 sinh 7 7ru; / 2 

At zero temperature, e(9) describes the excitation energy of the positively charged quasiparticles in the system which 
are the only ones filling the Fermi sea at V 7^ 0, T = 0. Taking a derivative of (^6|) we get the density of these 
quasiparticles 

p{6) J A d6> p{9'me e') = ±-e e . (29) 

The current then follows from a Boltzmann equation describing the one by one scattering of these charged excitation 
on the boundary. There are p(9)d9 such particles in [9, 9 + d9] and the probability that a particle of type a is reflected 
as one of type a' is (9 — log(Te))| 2 . This leads, following the arguments of fl to the current 

I(V,T B ) = e f A d9 P {9)\Rl{9 - log(T B ))| 2 . (30) 

J — 00 

Thus if we manage to determine A and solve the linear integral equation for p we get an exact evaluation for I(V, Tb). 
Since the integral equations are linear, this is possible by Wiener-Hopf techniques. 

In the previous equations, A is the Fermi rapidity defined by e(A) — 0; it is the edge of the Fermi sea. A 
can be determined, and the foregoing equations solved, by using a standard Wiener Hopf analysis. If we write 
K{6) = 6(9) - $(0), we have 



G 



with 



^)- 27r V 7 r(iW2)r(i/2 + iu,/2) 2 ' (32) 

where A is chosen to preserve the analyticity in the lower half plane. Here the conventions for the Fourier transforms 
are 



du) 



9(0) = J —e-^~g(u), ~g(cj) = J d0e^g{6). (33) 

Then following exactly the same steps as in || we find the solution for the fourier transform of the density 

i M+iui)A 
2h (1 + luj) 

The excitation energy of the particles follows 

e{u))e = -i— + — ; . (35) 

I OJ l OJ — I 

The condition e(A) = is equivalent to 

lim uji(ioj)e- lujA = 0, (36) 
which in turns gives the explicit value of the Fermi rapidity, A 

We are then left with the evaluation of the current given in ( p0| ) . After a few manipulations we get the exact expression 
This can be expanded in powers of V/Tg to get an IR expansion given by 

= 2 ^ Vf-1)" +1 T{2n/9p) n[2 ~ 9p) ' 3P 

h h (»-i)ir(i + =^)»[i + n(^)]V^ 

where we have reinstated the dependence on g p in the last expression and we have introduced the parameter 

T'b = Tb^- (40) 

for convenience. We can also separate the integral over rapidities to get an expansion in the UV. After a few 
manipulations we get 

«y,T B ) = ±NH)e*PP [ { V(-l)" +1 »" + 1 "r 1 ' 2 ,] (41) 

V ' ; 2h V ' J 2tt (1 + ioj) \^ Q K ' 2n/j + iuj 2 sin(i7rw 7 /2) J V ' 

with b = (Tee - " 4 ) 2 ' ^ 9p . We integrate the first sum by closing the contour in the lower half plane and the second term 
by closing in the upper half plane. In the second term, the zeroes of N(oj) cancels the poles of sinh(7TW7/2) and the 
one at oj = i, only the poles oj — 2ni/( r y + 2) with n integer contribute. We get 
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„2t/ % / rpi \ Ti(l-g p /2) 

e * TT I B 



I(V,T B ) = — g p + ^U n l-± , (42) 

n— 1 ^ ' 

with the coefficients 

U = l^V,_{-ir_ r(n gp /2) 

where again we have reexpressed things in terms of g p . From these two expressions we can check that there is an 
exact duality between the UV expansion and the IR expansion under the exchange g p — > 4/<7 p . This duality is more 
transparent through the relation 

I(V,T' B ,g p ) = ^-g p - ^I(V,T B A/g p ). (44) 

In order to make contact with perturbation theory, we need to find the relation between the boundary scale Tg and 
A. This is done in appendix B using a first order Keldysh computation. 



B. Conductance at finite temperature 



To compute the conductance in the case T^O and 7 an integer is difficult. The reason is, that the bulk scattering 
is not diagonal in the soliton antisoliton basis, and the TBA equations are written, in fact, for quasiparticles with 
no definite charge. The transport equation of jl2| becomes then ambiguous to use. In the case without spin, it 
was possible to conjecture a formula, based on limiting cases, that reproduced well numerical results, and later was 
established using more complicated functional relations. A similar conjecture in the present case would be 



d6 



2 J cosh (0 — In Tg/T) 



1 1 

1 _|_ e e 7 +l(0) 1 _|_ e e T +i(oo 



(45) 



It does go to zero in the infrared, as T/Tb — ► 0. In the ultraviolet, when T/Tb — > 00, it is easy to see that G goes 
to g p as desired. This is because, in that limit, the only contribution to the integral comes from the region 9 near 



infinity, where the e's go to a constant, in particular, e 



n(-oo) 



7+1 (while e 



1(00) 



7-1). 



The case 7 = 1/n, n an integer, belongs to the attractive regime of the sine-Gordon part of our scattering, and 
corresponds to the case with purely diagonal scattering. It is thus more favorable to study the conductance via the 
approach of E ]. However, while the attractive regime of the ordinary sine-Gordon model has been much studied, 
we are not aware of many such studies for the N = 1 sine-Gordon model of interest here. The complete S-matrix 
was found in |nj . The TBA has never been written, as far as we know. We will present a detailed discussion of this 
problem elsewhere. Here, we content ourselves by giving the relevant equations and the corresponding conductance. 
The diagram looks as follows 



' 6- 




n+l 



n-1 



FIG. 3. TBA diagram at finite temperature. 

The TBA equations have the form 

ej(e, V/T) =K*Y^ N ok In [l + e-^ k e 



lx kf> e k {9,V/T) 



(46) 



with K{9) — 2,i- coshne an( ^ J — 3, . . . , ri + 2, ^ n +i = ~ = w > Mfc = otherwise. In addition we have 

e a (0,V/T) =K* (In [1 + e £3 ] - In [l + e^ 1 ] ) 



(47) 
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and 



Here we have 



ei = -K* In [1 + e £2 ] -K'*]n[l + e~ £1 ] + 



sin 7r/2n 



K 



-, K> 



T cos7r/2n' 
cosh(n — l)iruj/2n 



2 cosh 7rw/ 2n ' 2 cosh 7rw/2n cosh 7rw / 2 

The other mass terms appear in the TBA in the form of asymptotic conditions: 



m 



(48) 



(49) 



(50) 



where nij = 2 sin ^'"^ j > 1, m„ + i = m„ +2 = 1. 

We can present, as a quick check of the validity of this TBA, a computation of the central charge of the theory. 
Setting x = e~^~°°^ we have the system 



1 =n(i+- 

Xj \ X\z 



N ]k /2 



supplemented by 



and 



1 

X 2 



Xi 



1 + 



x 3 



1/2 



(l + a;i)- 



-1/2 



1/2 

1+-) (1+Xl) 1 / 2 



The solution of this system of equations is 



1 



xi = 3; — = [j - - - 1, j = 2, . . . , n; 



1 1 

X„+i x n+2 2 



(51) 



(52) 



(53) 



(54) 



Similarly, introducing y = e e (°°) we clearly have that all the y's vanish, except y\ = 1. Therefore, using the general 
formula 



G 



where i is the dilogarithm function, wc find 



L 



y 



l + y 



2L 



k=l 



L(3/4) - L(l/2) 



where wc used the identities 



and 



LI 1 - 



1 + x 



L(x) + L(l -x) = 



The formula for the linear conductance follows 

G 



n r 

2J 



cosh 2 n{6 - \hTb/T) 1 + e e -+i( e ) 



(55) 



(56) 



(57) 



(58) 



(59) 



with the UV value G = 2 2 1 = .g p . In the limit Tb/T — > oo, the integral is dominated by the region — > oo where 
e„ diverges. It is then legitimate to simply expand the cosh in the integral, and one finds that G expands in powers 



of (ts - ) = {te) 3 " ' as rec l uircd - 
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IV. CONCLUSIONS 



To complete the analysis of this problem, we finally compute the boundary entropy. This is easy to do in the 
atractive regime. Since the boundary acts trivially on the kink degrees of freedom, the boundary scattering of the 
bound states is completely given in terms of the boundary scattering of the sine-Gordon part. This problem has been 
studied in (l^Jl^], and the R-matrices of breathers computed. The boundary free energy then follows 

Ti-4~2 n J/J 

3=2 J 7F 

where complete expressions for the Kj can be found in Jl6| (in the latter reference, A = n + 1). We can thus write 

n 

suv - s IR = lU) ln (! + Xj) + ( /(+) + /M ) Ml + (61) 

J=2 



where = J ^Kj{9). One finds thus 

n ; _ j 
SUV ~ SIR = ^ m 

J'=2 



(j ~ 1/2) 5 



C7 + 1/2)0" - 3/2) 



n in( n+1/2 )=-Hl + 2n). (62) 



2 \n-l/2j 2 



We find therefore suv — sjr = In y g 4 g > i n agreement with the IR fixed point being made of disconnected leads, and 

a computation of the boundary entropies using conformal partition functions |34j . 

In the repulsive regime, the computation is more difficult, again because the scattering on the impurity is non 
diagonal. However, when 7 is an integer, we expect, by analogy with the ordinary boundary sine-Gordon case, 

f*» = - T J £ cosh(.-ln(TB/T) ln ( 1 + e '" +1 ) (63) 
Hence we find the difference of entropies 

suv- sir = -ln(l + e- e ^- 00 A --hx(l + e - e ^°°A = - In-^- =bxj--- (64) 
2V / 2 V / 2 7 y g p g a 

as required. 

The tunneling problem with spin is expected to present many interesting features for arbitrary g p and g a . While 
it presumably is not always integrable, we have identified several integrable varieties besides the isotropic one just 
discussed. These include the case g„ + g a — 2, g p + q a = 4, and — + — = 2. We hope to report on the corresponding 
solutions (and, in some cases, new IR fixed points) in a subsequent publication. 
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APPENDIX A: ANALOGY WITH THE TWO CHANNEL KONDO PROBLEM, AND DERIVATION OF 

THE BOUNDARY S-MATRIX 

In this appendix, we discuss along more traditional lines how the impurity problem with g a = 2 can be mapped 
onto an anisotropic Kondo-type model. We start with a general anisotropic two-channel Kondo-type problem with 
hamiltonian 

/OO Q 
„, -°° ' •'' 
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where we have reformulated the problem so that all quantities are right movers. Here a are the spin indices and i 
indicates the "flavor" or fermion type. The impurity part is 

2 

Hi= Y, E J ^ A a>+(0)^(0) (A2) 

a,/3,i— 1 X—x,y,z 

where the a are Pauli matrices, while the S are generators of the U q sl{2) algebra 

[s+,s~] = <i 2S *-v- 2S ° (A3) 

q-q 1 

[S^S*] = ±S ± . (A4) 

The coupling constant are J x = J y ^ J z ; the parameter q is related with the anisotropy in a manner to be described 
below. 

Abelian bosonization (as in p3] which we follow closely here) is readily accomplished by introducing chiral bosons 

as 

i> ia (x) = == — . (A5) 

x/ZTTCL 



Jz ^ R 



In dx 



The bulk part becomes just a free right moving boson hamiltonian, while the impurity part reads 

+ L ( 5 + e iVSF^(o) +S - e - iV ^^}co8V^^ f (0). (A6) 
na I > 1 

In (A6) we have expressed the result in terms of new boson fields defined by 

€ = \ E & 

1 

<£ = \ E <«<t>?a (A7) 



2 

i,a— 1 



2 

Z -6 R 

it uv 



l.OL — 1 



1 2 

2 E VaaVii^ia 



for the charge, spin, flavor and spin-flavor fields respectively. We find it convenient to absorb the J z term by performing 
a unitary transformation with U = exp[i-^S z (j) R (0)] that suppresses the J z term in the impurity hamiltonian which 
reads now 



Hl = + ± ( 5 + eiv ^^(o)/V2 + s - e -iy^s?^lo)/V5\ cos ^« (o). (A8) 



.]_ 

na 

with = ^/n — . Normalizations are such that the dimension of the impurity operator is d = i + Sr. . I n 

what follows, we shall use g p (instead of J 2 ) as the only dimensionless parameter of the problem. The quantum group 
parameter for the impurity spin is then q — e md . 

The boson fields <j) c and <pf are totally decoupled from the hamiltonian. If q is a root of unity, and if we chose for 
the impurity spin a periodic (cyclic) representation of the quantum algebra, then, by following the argument of [fl2f , 
it is easy to see that the problem is equivalent to the "double" impurity sine-Gordon model 

Hl = — cos(^4^</>f(0)/V2)cosV4^/)f f (0). (A9) 
na 1 

with J oc J' . Equivalently, we can fold the problem to recover left and right movers on a half line, transforming in 
this way an impurity into a boundary problem 
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J 



Hi = — cos(v/2^~^ s (Q))cos2V7r^ s/ (0). (AlO) 
ira 

This coincides with the hamiltonian in the main text, after identification P = 4> s and 9 a = <j) s f. 

On the other hand, the integrability of the higher spin, multiflavour Kondo hamiltonian is well established. Al- 
though the isotropic case is usually considered, the proof immediately generalizes to an impurity spin in an arbitrary 
quantum group representation [p5|-p7f. To understand the integrability structure, it is useful to perform non abelian 
bosonization. Introduce the currents 

J x = ^tA^a (AH) 

for the U(l) charge, SU(2) spin and SU(2) flavor currents respectively. The bulk hamiltonian is then quadratic in 
currents 

H a = JJ + J x J x + J a J a . (A12) 

In the same manner the interaction with the impurity reads 

Hi = J(J+S-+ J~S+) + J Z J Z S Z . (A13) 

Only the spin currents are interacting with the impurity, and we can forget in what follows the free charge boson (a 
c = 1 theory) , as well as the flavor current (a c = | theory). 

The massless scattering description of the general Kondo model was given in [ p8| . The bulk degrees of freedom 
are the massless limit of the well known N = 1 supersymmetric sine-Gordon model. The effect of the impurity is 
described by an R-matrix, whose form depends on the impurity spin s. In the underscreened case (s > 3/2), it is 
given by a solution of the Yang Baxter equation corresponding to scattering a spin 1/2 (the (it, d) doublet) through a 
spin s — 1 (this renormalisation of the spin occurs, because the scattering theory is realy an IR description, and that, 
in the IR, 2 electrons screen the impurity). The RSOS degrees of freedom scatter then trivially. On the other hand, 
the cyclic representations behave in many ways like an infinite spin. As a result, the R-matrix is given by an object 
similar to the underscreened case, with the (m, d) doublet scattering through a cyclic representation of parameter 
q = — e~™/ 7 . As discussed in Jl^ ], this cyclic spin can be gauged away, and one obtains the well known boundary 
sine-Gordon R-matrix for the (u, d) degrees of freedom, as discussed in the text. 



APPENDIX B: KELDYSH COMPUTATION. 



In order to compute the differential conductance we need to use the Keldysh formalism since the system is driven by 
reservoirs. To do so, we will use the formulation of the model on the full line as described by the equations (|Io|). The 
effect of the reservoirs can be implemented by shifting the charge fields, 8 P — > O p + \J g p /na{t)/2 with dtait) = V. 
Under this prescription the current is evaluated by taking the functional derivative 

1(t) . -.MM, (Bl) 

where we have used conventions in which h = e = 1. Using the Keldysh contour, C, which goes from — oo to co and 
then comes back (see figure 3), to expand the partition function we obtain to first order 



oo 



-oo C 

FIG. 4. Contour for the perturbative evaluation. 
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7 (2) (0) = _i^iP jf A sin[|a(i)] {P'X^+Pt^A-t)} (B2) 

where /i = ± depending on the location of t, ie upper part of the contour or lower part of the contour. The functions 
P^ 41 (t) are the corresponding contraction of the vertex operators time ordered on the Keldysh contour 

p ^ = Js±W^> p * ±T W = (*iW (B3) 

To this order the result can be found explicitely, it is given by 

J "'°» = - 2r[ tep + i ,,)/2] (f) 5 <B4 > 

This is in agreement with the Bethe ansatz solution given in the bulk of the text provided we make the identification 
(putting h — 2-7T and e = 1 in the TBA expressions and g a = 2 in the previous perturbative expression) 

r(gp/2 + i) V (2~g P \ 9p/2 T i B - 9P /^ (B5) 



5p(2-5p) V r (5p/4) / V n 9 P 
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